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SU(3) Flavor Breaking in Hadronic Matrix Elements for B — B 
Oscillations 
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^ Department of Physics, Washington University, St. Louis, MO 63130 
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We present an analysis, using quenched configurations at —5.7, 5.85, 6.0, and 6.3 of the matrix element 
Mhi = {Phi\h'yi_i{l~"f5)lh'yfj,{l~^5)l\Phi) for heavy-light pseudoscalar mesons. The results are extrapolated to the 
physical B meson states, B'^ and _Bg. We directly compute the ratio Mbs/M^d, and obtain the preliminary result 
Mbs/Mbd ~ 1.54(13)(32). A precise value of this SU(3) breaking ratio is important for determining Vm once the 
mixing parameter Xs for B^ — B^ is measured experimentally. We also determine values for the corresponding B 
parameters, i?i,s(2GeV) = Bbd(2GeV) — 1.02(13), which we cannot distinguish in the present analysis. 
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1. Introduction 

We present a direct calculation of the AF ~ 2 
heavy- light mixing matrix element, 

Mw(m) = {Philhpa - 75)^^(1 - 75)^1^^/./), (1) 

where we have suppressed the scale dependence of 
the local four quark operator. As is well known, 
these matrix elements govern B^—B^ and B^—B^^ 
oscillations |jl|,||. In the above h and I denote 
heavy and light quark fields, P^i the correspond- 
ing pseudoscalar meson, and fi is the energy scale 
appropriate to the calculation. In particular, we 
present a first direct calculation of the SU(3) fia- 
vor breaking ratio. 



Tsd = Mbsip)/Mbdip) 



(2) 



Our preliminary result is — 1.54(13)(32) 
where the first error is statistical and the sec- 
ond is from uncertainty in extrapolating to the 
B mass and to a — > 0. The importance of this ra- 
tio is that, in conjunction with the experimental 
measurement of — B'^ oscillation (when that 
becomes available), it should allow the cleanest 
extraction of the crucial CKM parameter Vtd- We 
note that the above value of is lower than the 
result reported at the conference (1.81) which is 
due primarily to a change in how we extrapolate 
to the B: the difference gives a systematic error. 



Presently, Vtd is deduced from B'^ — oscilla- 
tion vialHl 



Xbd 



(AM) 



bd 



oc mldBbditJ-)fbd\ytd\ 



(3) 



where rubd, L^ = Tbd, and fbd are the mass, 
life time, and decay constant of the B'^ meson, 
and {AM)bd is the mass difference of the two 
mass eigenstates of the i?° — i?" system. Xbd 
is the mixing parameter characterizing the oscil- 
lation and has been determined experimentally, 
Xbd = 0.71(6) Bbd is the so called bag, or B, 
parameter. To extract Vtd from Eq.( ||) requires 
knowledge of two hadronic matrix elements, fbd 
and Bbd- These are being calculated by using 
lattice and other methods, fbd may eventually 
even be measured experimentally through, for ex- 
ample, the decay B TVr- However, Bbd is a 
purely theoretical construct which is inaccessible 
to experiment. Thus determination of Vtd from 
experiment through use of Eq. ^ will ultimately 
be limited by the precision of the nonperturbative 
quantity fbd^bd- These parameters are related to 
the matrix element, Eq. (^, via 



Mbdin) = -fbdmldBbd- 



(4) 



Now making the replacement — > s in Eq. (|^) 
and taking the ratio with Eq. (0), we arrive at an 
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alternate way to extract Vtd, 



\Vtd\ 



Tsd- 



Tbs Xbd 
Ud Xbs 



(5) 



Note that Vts in Eq. (||) is related by three genera- 
tion unitarity to Vcb and is therefore already quite 
well determined, \Vts\ w \Vcb\ = 0.041 ± 0.003[|]. 
The important distinction between using Eq. (H) 
instead of Eq. (||) is that the former requires only 
knowledge of corrections to SU(3) flavor symme- 
try while the latter requires the absolute value of 
the matrix element Mbd- It is also important to 
realize that since rsd is a ratio of two very similar 
hadronic matrix elements, it is less susceptible to 
common systematic errors in lattice calculations, 
including scale dependence, matching of contin- 
uum and lattice operators, and heavy quark mass 
dependence. Indeed, the ratio r^d is, to an excel- 
lent approximation, renormalization group invari- 
ant, even though the individual matrix elements 
Mbs and Mbd are scale dependent. 

2. Simulations and Results 

Table ^ summarizes our quenched lattices and 
the valence Wilson quark hopping parameters, ni 
and Kh, used to construct quark propagators. For 
each Ki and Kh in Table |l| we calculate a quark 
propagator using a single point source at the cen- 
ter of the lattice and a point sink. 

In Fig. 0we show example results at6/g^ = 6.3 
for Mhi versus Kh for each value of ki. The lattice 
matrix elements are found through simultaneous 
fits to the three point function corresponding to 
Eq. |l| and the two point function of the corre- 
sponding heavy-light interpolating operator. In 
Fig. |l] and the following, we have already matched 
the lattice operator to the continuum one by using 
the one loop perturbative result from Refs. [^j?). 
The scale dependent renormalization Z factor is 
calculated in the NDR scheme at scale ^ — 2.0 
GeV. The coupling is tadpole improved and eval- 
uated at scale 1/a, and we include the KLM nor- 
malization for heavy quarks. 

Results for the physical B and Bg meson sys- 
tems follow from a series of fits to the lattice data 
which we use to extrapolate in the two parameters 
Kh and Ki. Since the data are highly correlated. 



we use covariant fits and a jackknife procedure 
at each step to account for the correlations. We 
take the form of our fits from chiral perturbation 
theory and expectations based on heavy quark 
effective theory (HQET). 
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Figure 1. The matrix element Mhi at — 6.3. 

Kci&lg^) and KsiQ/g"^) are determined from ei- 
ther linear or quadratic fits to rnfi, as a function of 
K^^ and KjT^ {I' refers to the strange quark). The 
values for Kc and Kg are summarized in Table |^. 
Finding Ks requires the scale a, which we set from 
o-Itt, to determine the lattice value of the kaon 
mass arriK (a-i is tabulated in Table (|)). We 
note that at = 5.7 the choice of the coupling 
constant scale for Za, the axial current renor- 
malization, has a significant effect on the lattice 
spacing determination; Za differs by ~ 7% when 
the scale changes from 1/a to n/a. 

Next, using chiral perturbation theory for 
heavy- light mesons |§,|l^, we extrapolate Mhi to 
Ki = Kc- We do not include chiral logarithms 
in our fits since the light quark masses used in 
the extrapolations are relatively heavy, k/ « k^. 
The results for Mm at 6/5^ = 6.3 (see Fig. |) 
show a smooth linear behavior. The confidence 
levels for the extrapolations are much lower for 
Q/g^ = 6.0(24'^) than for the other points; in 
addition, the effective values of Mhi show poor 
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Table 1 



Summary of simulation parameters. 




conf. 


size 






5.7 


83 


16^ X 33 


0.160 0.164 0.166 


0.115 0.125 0.135 0.145 


5.85 


100 


20^ X 61 


0.158 0.159 0.160 


0.107 0.122 0.130 0.138 0.143 


6.0 


60 


16^ X 39 


0.152 0.154 0.155 


0.103 0.118 0.130 0.135 0.142 


6.0 


100 


X 39 


0.152 0.154 0.155 


0.103 0.118 0.130 0.135 0.142 


6.3 


60 


243 X 61 


0.148 0.149 0.150 0.1507 


0.100 0.110 0.125 0.133 0.140 



plateaus, so we exclude this point in our final de- 
termination of r^d- 
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Figure 2. M^i at 6/(7^ = 6.3 extrapolated to 
Ki = Kc- Curves for the other values of Q/g^ are 
similar. 

Inspired by HQET, we continue by fitting Mhc 
to a polynomial in the inverse heavy meson mass, 
m'^l , and then extrapolating to the B meson 
mass (see Fig. ||). Again, we use /^r to set the 
scale. Of course, the same procedure can be car- 
ried through for the heavy-strange mesons by first 
extrapolating the data to Kg instead of Hc- The 
resulting curve is also shown in Fig. |^. Fits which 
include all mass values (dashed line) generally 
have low confidence levels. Moreover, the lighter 
points have smaller statistical errors and domi- 
nate the fits, yet they are far from the B mass. 
We therefore omit the lightest two points at each 
value of 6/(7^ and use a completely constrained fit 



Table 2 



Inverse lattice spacing and critical and strange 
hopping parameters. 





a-^(GeV) 


tie 




5.7 


1.45(10) 


0.16969(10) 


0.1642(7) 


5.85 


1.64(20) 


0.16163(9) 


0.1577(9) 


6.0 


2.06(17) 


0.15711(7) 


0.1548(4) 


6.0 


2.08(13) 


0.15714(4) 


0.1543(4) 


6.3 


3.37(47) 


0.15226(16) 


0.1506(4) 



to extrapolate to the B. Compared with fitting 
all points, this results in systematically lower val- 
ues of rsd, changing the central value from 1.81 
(which was reported at the conference) to 1.54. 
We take the 0.27 difference as the systematic er- 
ror of the extrapolation. The ratio M^s/Mbd is 
shown as a function of a{Q/g^) in Fig. |^. At 
6/(7^ = 5.7, our heaviest mass is still quite far 
from the B mass, so we also ignore this point in 
the extrapolation to the continuum limit. 

Generally, we expect the Wilson quark action 
to introduce discretization errors of order a in all 
observables. However, in a ratio of two similar 
quantities, we might expect a large cancellation 
of the lowest order discretization errors. A con- 
stant fit gives Mbs/Mbd — 1.54(13) while a linear 
fit gives 1.72(67). Since the coefRcient of the lin- 
ear term is only 0.3 sigma from 0, we quote the 
constant fit result as our central value and use 
the difference as an estimate of the systematic 
error of the continuum extrapolation. Adding 
that error in quadrature with the systematic er- 
ror from the extrapolation to the B mass, we get 
Tsd = 1.54(13)(32). 

The extraction of the individual values of Mbd 
and Mbs is clearly expected to have larger errors. 
Conventionally (|Jl|,|| these matrix elements are 
given in terms of the corresponding B parameter 
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Figure 3. Mhc (octagons) and Mhs (squares) as a 
function of the inverse heavy-down(strange) me- 
son mass, at Q/g^ = 6.0. The dashed line shows 
the effect of the Ughtest points on the fit. 

defined in Eq. (^) . Carrying out the above fitting 
procedure for B^,d{^^), we find a constant fit yields 
5^(2 GeV) — 0.97(3) while linear extrapolation 
gives 1.02(13). We cannot however distinguish 
Sbs(2GeV) from i?M(2GeV) since our data for 
Bhi versus are fit equally well to constant or 
linear forms. This was not true for Mhi as is 
evident from Fig. ^ Taking the linear result, we 
quote BM(2GeV) = Bbs{2GeV) = 1.02(13). 

We recall that until now [0J§], lattice results 
for the SU(3) breaking ratio rsd have been ob- 
tained by using Eqs. (||) and (|^) and the lattice 
measurements of fbd{s) and Bi,d(s)- A reasonable 
estimate is hs/hd « 1.13 ± .10 [g||,|ll[ (we 
are presently calculating this ratio on our lat- 
tices). As indicated above, the ratio of B pa- 
rameters is consistent with unity, and the ratio of 
masses is 1.017 jj]. Therefore, the conventional 
method leads to Tgd ~ 1-32 ± .23 compared to 
1.54 ± .13 ± .32 obtained with our direct method. 
Thus the two methods are quite consistent. How- 
ever, as we have emphasized, the direct method 
offers many distinct advantages, and future lat- 
tice computations should be able to improve the 
precision in our determination of the ratio r^d- 




Figure 4. The SU(3) flavor breaking ratio 
Mbs/Mhd versus the lattice spacing a. The points 
denoted by crosses were used in the fit (solid line). 
The circles (/? = 5.7 and (3 — 6.0, 24'^) were omit- 
ted, for reasons explained in the text. The burst 
shows Tsd extrapolated to a = 0. 
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